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Abstract: We suggest a real right-handed sneutrino, N\, as a good dark matter candidate 
in a supersymmetric Z' model realizing the seesaw mechanism. When the extra gaugino, 
Z' , is lighter than Z' , the thermal freeze-out of the dark matter annihilation to right- 
handed neutrinos, N±Ni — > NN, through the t-channel Z' exchange is shown to produce 
the right dark matter density. It is essential to include the decay and inverse decay of N 
in this process, otherwise N decouples too early and thus dark matter is overproduced. At 
the LHC, the search for the seesaw mechanism can be made by observing the signatures 
of pp — > Z'Z' — > NN+ p"t as Z' can be copiously produced from the cascade decays of 
gluinos/squarks, which is complementary to the search of pp — >■ Z 1 — >• NN. This may also 
open up a promising new channel of finding the Higgs boson from the displaced N decay. 
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1. Introduction 

Supersymmetry at the TeV scale has a nice feature that it provides a good dark matter 
(DM) candidate, the lightest supersymmetric particle (LSP) which is stable when R-parity 
is imposed for the proton stability. A typical LSP in the supersymmetric Standard Model 
is a neutral gaugino with a small mixture of Higgsinos whose relic density is naturally in 
the right range due to its weak interaction property [1]. The supersymmetric Standard 
Model may be extended to include a new sector in which the LSP resides. The observed 
neutrino masses and mixing require such an extension in various ways [2]. A standard 
example is to introduce right-handed neutrinos which are singlets of the Standard Model 
(SM) gauge group and have only Yukawa interactions with active neutrinos. Then, a right- 
handed sneutrino, which is the scalar superpartner of a right-handed neutrino, can be the 
LSP and thus a dark matter candidate. However, the cosmic production mechanism of the 
right-handed sneutrino dark matter can be different from the usual neutralino. 

If neutrinos are Dirac particles without invoking the seesaw mechanism, the required 
neutrino Yukawa coupling is of order 10 -13 , and thus the right-handed sneutrino can never 
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be thermalized. However, thermal regeneration through the tiny Dirac neutrino Yukawa 
coupling can be effective to produce a sizable abundance of the right-handed sneutrino [3] . 
If the seesaw mechanism is realized at the TeV scale, one has the neutrino Yukawa coupling 
of order 10 -7 which, then, leads to the overproduction of the right-handed sneutrino LSP. 
To avoid this, there must be a thermalization process through which its relic density can 
be suppressed appropriately. One of the earliest proposals for such a thermalization is 
to assume an unconventionally large A term leading to a sizable mixing between the left- 
handed and right-handed sneutrinos [4]. More recently, the non- minimal supersymmetric 
Standard Model was considered to allow a coupling of the right-handed neutrino with 
an extra singlet field introduced in the Higgs sector [5]. Note also that one can invoke 
the inverse seesaw mechanism which allows for a possibility of order-one neutrino Yukawa 
couplings [6]. 

In this paper, we propose a new way of generating thermal freeze-out density of the 
right-handed sneutrino dark matter in U(l)' seesaw models which may arise from a grand 
unification group [7]. In the case of Dirac neutrino models with no lepton number viola- 
tion, the right thermal abundance was shown to be generated through the Z' resonance 
enhancement of the right-handed sneutrino annihilation [8]. However, this feature is not 
shared with U(l)' models realizing the seesaw mechanism in which (complex) right-handed 
sneutrino components are split into two real scalar fields, an LSP and a heavier state, with 
a large mass difference driven by the lepton number violating soft (Majorana) mass term 
and they couple non-diagonally to the Z' gauge boson. Thus, the Z' interaction is typically 
too weak to prevent the overproduction of the right-handed sneutrino LSP. 

It will be shown that the right-handed sneutrino annihilation to right-handed neutrinos 
through the i-channel U(l)' gaugino, Z', exchange can be efficient enough to produce the 
observed dark matter density. In this process, the decays and inverse decays of the right- 
handed neutrinos through the small Yukawa coupling play an important role in keeping 
the right-handed neutrinos in thermal equilibrium and thus controlling the dark matter 
relic density, which is a distinguishable feature of the thermal history of the right-handed 
sneutrino dark matter compared with the conventional neutralino LSP dark matter. From 
the numerical analysis of the Boltzmann equations, we will show how the dark matter 
freeze-out density arises depending particularly on the Z' mass and the effective neutrino 
mass corresponding to the dark matter sneutrino Yukawa coupling. 

Searching for an extra gauge boson Z' through dilepton resonances is one of the key is- 
sues in collider experiments [9]. A novel feature of U(l)' seesaw models is that right-handed 
neutrinos can be produced directly from the Z' decay and their subsequent decay leaves a 
clean signature of same-sign dileptons appearing at sizable displaced vertices depending on 
the corresponding neutrino mass [10, 11]. Observing these signals at the LHC will indicate 
the Majorana nature of neutrinos confirming the seesaw mechanism [12]. Furthermore, 
it provides an interesting possibility to discover the Higgs boson through displaced 6-jets 
[13] coming from the right-handed neutrino decay. The supersymmetric model with the 
right-handed sneutrino LSP enjoys these features with the production of Z' and its decay 
to the right-handed neutrino and sneutrino dark matter. Independently of the Z' produc- 
tion, Z' can be produced copiously from the cascade decays of gluinos and squarks, and 
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thus the search for the presence of the Majorana right-handed neutrino can be made more 
effectively. In particular, it is expected to have a large number of events for the same-sign 
dilepton or four displaced 6-jets plus large missing energy from pair-produced Z' if it is the 
next LSP. 

This paper is organized as follows. In Section 2, we introduce a U(l)' seesaw model, 
and analyze the masses of Z' and Z' and interaction vertices relevant for our calculation. 
In Section 3, the thermal abundance of right-handed sneutrino dark matter is calculated 
for various model parameters such as the masses of Z 1 and the right-handed (s)neutrino, 
the neutrino Yukawa coupling, and the U{1)' gauge coupling. In Section 4, LHC signatures 
of the model are discussed focusing on the pair production of the right-handed neutrinos 
from the decay of Z' and Z 1 , which can provide a test for the seesaw mechanism and a new 
channel for the Higgs production. We conclude in Section 5. 

2. A supersymmetric U(l)' seesaw model 

Among various possibilities of an extra gauge symmetry U(l)' and the presence of the 
associated right-handed neutrinos [7], we will take the U(l) x model for our explicit analysis. 
The particle content of our U(l) x model is as follows: 



5*7(5) 


10 F 5 F 1(JV) 


5h 5h 


l(X) 


l(5i) 1(5 2 ) 


2y/WQ' 


-1 3 -5 


2 -2 





10 -10 



(2.1) 



where SU{5) representations and U (1)' charges of the SM fermions (10^, 5f), Higgs bosons 
(5#,5#), and additional singlet fields (N, 5, S' 12 ) are shown. Here N denotes the right- 
handed neutrino, X is an additional singlet field fit into the 27 representation of Eq, and 
we introduced more singlets 5^2, vector-like under U{\)' , to break U(l)' and generate the 
Majorana mass term of N. Note that the right-handed neutrinos carry the largest charge 
under ?7(l)y an d thus the corresponding Z' decays dominantly to right-handed neutrinos. 
Furthermore, the additional singlet field X is neutral under U(l) x so that it can be used 
to generate a mass to the U{\)' Higgsinos as will be discussed below. The gauge invariant 
superpotential in the seesaw sector is given by 

W seesaw = VijUHuNj + ^SiNiN, , (2.2) 

where Li and H u denote the lepton and Higgs doublet superfields, respectively. After the 
f7(l)' breaking by the vacuum expectation value {Si), the right-handed neutrinos obtain 
the mass = AjVj(5i) and induce the seesaw mass for the light neutrinos: 

m-j = -VikVjk-^ ■ (2.3) 

To generate the atmospheric neutrino mass fh u = 0.05 eV with the right-handed neutrino 
mass scale of ~ 100 GeV, we need the Yukawa coupling of y v ~ 4 x 10~ 7 . For the 
breaking of U{1)' , one can consider the following schemes. The first option is to introduce 
the /j! term [14]: 

W' = ti'SiS 2 , (2.4) 
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and arrange radiative U(l)' breaking with a large Yukawa coupling, say Atv 3 , mimicking the 
electroweak breaking in the Minimal Supersymmetric Standard Model (MSSM). Another 
possibility is to consider 

W = \XS 1 S 2 + ^X 3 (2.5) 
as in the non- minimal Higgs sector [15]. The corresponding scalar potential is given by 



V = V susy + Fsoft + V D 



where 



susy 



Kioft 



V D 



=X,S 1 ,S 2 ,N l ,H u ,L l 

y v A v LH u N + ^A N SiNN + XA s XSiS 2 + ^A K X 3 + h.c. 



+m 2 x \X\ 2 + m 2 s [\Si\ 2 + |5 2 | 2 ] + m 2 ~\N\ 2 , 



J2 



Q' Sl \S 1 \ 2 + Q' S2 \S 2 \ 2 + Q' N \N\ 2 + 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



In this case, the /jf term is generated through the vacuum expectation value of X: fjf = \{X) 
which, in the leading term, is given by 



s I 



(2.10) 



where A K is the trilinear soft term and ms is the soft mass of 5" [15]. The vacuum expec- 
tation values of S\ t 2 can be developed radiatively for A ~ 1. 

Given the vacuum expectation values (Si^), one gets the Z' gauge boson mass: Mg, = 
2g' 2 ^2iQ's.{Si) 2 where g' denotes the U(l)' gauge coupling. At the moment, the most 
stringent bound on the Z' mass comes from the electroweak precision data. For the U{\) x 
model, we get Mz> > 1.14 TeV for the reference gauge coupling g' = ^Jhj 3g 2 tan 6\y ~ 0.46 
[16], which gives (Si) > 800 GeV for (Si) = (S 2 )- Note that the value of the gauge coupling 
g 1 can be scaled down depending on the particle content below the grand unification scale. 
For our numerical analysis, we will take the above reference value. 

Important mass parameters for our analysis are the U(l)' gaugino mass and the dark 
matter mass. Defining tan/3' = (S 2 )/(Si), the U(l)' gaugino-Higgsino mass matrix in the 
basis of [Z', Si, S 2 ] is given by 



M 



Mz>cp 











(2.11) 



where thm is the soft supersymmetry breaking mass. In the limit of // 3> mM,Mz', the 
lightest state has the mass m z , given by m z , ~ tyim + M z ,s 2 ^ij^. In the following, we 
will work in this limit. As a light Z' is favored for our dark matter relic density, we will 
work in this limit taking free parameter without bothering the specific values of 



- 4 - 



tan/3' and //'. Among three right-handed neutrino superfields, let us take the lightest one 
N (we will suppress the flavor index in the following discussions) whose scalar component 
contains the LSP. Denoting its supersymmetric mass by m^, the right-handed sneutrino 
N has the mass terms: 

Vmass = {m 2 N + m%- ±M%,C2fp)\N\ 2 - ^B N m N (NN + N*N*) , (2.12) 

where B^m^ = —\nAn{S\). Due to the lepton number violating (Majorana) mass term 
of Bntun which is assumed to be positive, the real and imaginary components of the 
sneutrino, N = {Ni+iN 2 )/^, get a mass splitting and their masses are given by = 
m 2 N + — =p B^niN- The real scalar field N± is taken to be the LSP dark matter. 

2.1 Couplings of N and N\ with Higgses 

From the superpotential in Eq. (2.2), simplified here as W seesaw = y v LH u N +^ m N N N , and 
the soft terms in Eq. (2.8), one can find the couplings of N and N involving leptons/sleptons 
and the Higgs bosons. For this, one uses the usual transformation to the Higgs mass 
eigenstates as follows: 

H% = -j={v d + h° d + ia§ , H° u = ^={v u + h u + ia° u ) , 

h d = cos aH — sin ah , h^ = sin aH + cos a h , 
a° d = sin /3 A - cos /3 G° , a° u = cos /3 A + sin /3 G° , 
H% = sin/3fl' ± -cos /3G ± , = cos /3 il ± + sin /3 G ± , (2.13) 

where tan/3 = v u /vd and tan 2a = tan 2/3( ^f^^f ). 

Then the right-handed neutrino ./V has the coupling with the left-handed neutrino v 
and charged lepton I: 

C NL = ^= Nv (h cos a + H sin a + iA cos /3) - y u NIH + cos f3 + h.c. (2.14) 
a/2 

in the Weyl fermion notation. Due to the vacuum expectation values of the Higgs bosons, 
there arises mixing between N and v and the corresponding mixing angle is given by 

e Nv « . (2.i5) 

This mixing induces the N—v—Z and A-Z-W -1- interaction leading to the A decay to the 
usual gauge bosons. For completeness, let us write down the scalar interaction of N = 
(N x + iN 2 )/V2: 

Vm = y u m N DH° u N* + y v u{A L H° u - »H° d *)N (2.16) 
-y v m N lH+N* - yJ(A L H+ + fiHf)N + h.c. , 

from which one can read off the couplings in terms of the mass eigenstates. The above 
scalar interaction leads also to mixing between A and v whose mixing angle is given by 

a _ (m N ± A L T V cot /3)v u 

°Ni**hj-V» ^2 (m 2 -ml ) ■ {IA<} 
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2.2 Decay of the right-handed neutrino N 

From the couplings (2.14) and the mixing (2.15), various channels open up as the de- 
cay modes of the right-handed neutrino, N. The corresponding decay widths are given as 
follows: 

r(AT -> uh) = r (N -> vh) = ( 1 - ^$ V cos 2 a , (2.18) 

8 8tt \ m N J 



iV 



8 8ir \ m 2 N 

9 / 2 \ 2 



r(iv^^+) = r(iv^^) = i^^fi-^ cos 2 /3 

4 8^ V ™?v / 



lylm N / M 2 \V , M 



r<* -+ „Z) = T(N -> m = - = 1 - ^ 1 + 2^ ) sin 2 , 



Ar/ \ '"-AT 
2 



r(AT -> IW+) = T(N -> W") = 1 ^ f 1 - ^ f 1 + 2^ sin 2 (3 . 

8 8vr V m N J \ m N J 

To quantify these decays, we will use the notion of the effective neutrino mass defined by 

rhu = |^ (2.19) 

which can be of order of the observed neutrino mass or smaller. 
2.3 Couplings of Z' and Z' 

The U(l)' gauge boson and gaugino couplings in our model is summarized in Table 1. 
For the Higgs-neutralino(chargino)~Z / coupling, the diagonalization matrices of the usual 
neutralinos and charginos are defined by 

H° u = N 4l \°, H2 = N 3i \1 (2.20) 
H+ = V 2i \f, Hj = U 2i \r 

in the two component notation. Moreover, the four component mass eigenstates are defined 
in terms of the two component fields as 



3. Relic density of the right-handed sneutrino dark matter 

In this section, we will calculate the thermal relic abundance of the right-handed sneutrino 
DM N\. Since the DM particle N\ cannot annihilate through the Z' interaction, one may 
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Table 1: Couplings of Z' and Z' in the U (l) x model where g' and Qx are the U(l) x gauge coupling 
and charge of the field X, respectively, and denotes the 4- momentum difference between two 
bosons in the vertex. 



consider its annihilation to right-handed neutrinos through the i-channel Z' exchange. 
In this case, the question is whether the right-handed neutrino iV remains sufficiently in 
thermal equilibrium. Consequently, it is very crucial to study the interaction rate for 
reactions of N. As will be shown in detail, the Z 1 interaction of N alone is too weak to 
give the right dark matter relic density but the inclusion of the decay and inverse decay 
of N can significantly change the result. For the following analyses, we will assume the 
numerical values of rrih = 115 GeV, tan/3 = 10, and g' = 0.46 unless otherwise stated. 



3.1 Freeze out of N\ 

The annihilation cross section of the DM candidate N\ through the process, N\N\ —> NN, 
can be large enough to suppress arbitrarily its relic density depending on mass parameters 
such as ,rriN, and m^, as far as N is in thermal equilibrium long enough. However, 
the interaction of N is either through the heavy gauge boson Z' or the tiny neutrino 
Yukawa coupling y v and thus might be too weak to keep N in thermal equilibrium during 
the process of freeze-out of N\. Thus, to study the thermal history of JVi through the 
annihilation in Figure 1(a), one has to consider also the evolution of N determined by its 
annihilation (Figure 1 (b)) and decay (Figure 1 (c), (d)). 

For the analysis of the relic abundance of the N\ , let us set up the following coupled 
Boltzmann equations for the evolution of the number density rii of particle i: 



dn N, 

dt 



-3Hn^ - {vfrVfij 



9m 

gN 



n N 



(3.1) 
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Figure 1: Annihilation and decay channels of N± and N. In panel (d), <j> = h, H and A. 



V 9N 



9n, 



= -3Hn N - {(Jnvn) [(n N ) 2 - 0^) 2 ] + 

- T N (n N - n^) , (3.2) 

where i? is the Hubble parameter, and n^ q , and gi are respectively the equilibrium 
number density, relative velocity and number of internal degrees of freedom of particle i. 
The first terms on the right-hand sides of Eqs. (3.1) and (3.2) account for the dilution due 
to the expansion of the Universe. The second term on the right-hand side of Eq. (3.1) 
describes the forward and backward reactions of N±N\ annihilation to NN through the 
4-channel Z' exchange (Figure 1 (a)). The second term on the right-hand side of Eq. (3.2) 
refers to the forward and backward reactions of NN annihilation to the SM fermion pairs 
// through the s-channel Z' exchange (Figure 1 (b)). The third term represents the 
effects of N\Ni annihilation to AW through the t-channel Z' exchange (Figure 1 (a)). The 
last term describes the decays and inverse decays of the right-handed neutrinos shown in 
Figure 1 (c), (d). 

In solving the Boltzmann equations (3.1) and (3.2), it is useful to introduce the variable 
Yi = m/s describing the actual number of particle i per comoving volume, where s is the 
entropy density of the Universe. Solving the coupled differential equations (3.1) and (3.2), 
one can find Yi as a function of x = m^jT . 

In order to see the effect of the decay more clearly, let us first present the solutions 
of the coupled Boltzmann equations excluding the last term on the right-hand side of 
Eq. (3.2) which accounts for the decays and inverse decays of the right-handed neutrino. 
The result is shown in Figure 2, where we use the following mass parameters: tu^ = 300 
GeV, m N = 260 GeV, m ? , = 600 GeV, and M z > = 1200 GeV. Since N x can remain in 
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Figure 2: The actual number of N± and N per comoving volume without the decay effect term. 
Blue solid and red dashed lines show = /s and Yjy = / s respectively. The other 
parameters are fixed as follows: = 300 GcV, m N — 260 GeV, mg, = 600 GeV, and M z * = 1200 
GeV. 

thermal bath through the interactions of N, the freeze out temperature of Ni is determined 
by that of N, if the interactions of N are weaker than those of Ni. Therefore, without the 
decay effect term, N\ is decoupled from thermal equilibrium when N is decoupled as can 
be seen from Figure 2. 

The right-handed neutrino mostly annihilates to SM fermion-antifermion pairs // 
through the s-channel Z' exchange, but the cross section for this process NN — > ff is sup- 
pressed by large Z' gauge boson mass, Mz' > 1.14 TeV [16]. Moreover, in the zero- velocity 
limit, this annihilation cross section of right-handed neutrinos is more suppressed since the 
annihilation cross section for this process has no s-wave component. Thus, the annihila- 
tion of the right-handed neutrino usually freezes out earlier than that of the right-handed 
sneutrino DM candidate. Consequently, the right-handed sneutrino N% is generically over- 
produced due to the early decoupling of N as shown in Figure 3. 

The decays and inverse decays of right-handed neutrinos through the Yukawa coupling 
y u turn out to play a crucial role in evading the previous overproduction problem. In 
Figure 4, we present the evolution of the number density of N\ and N per comoving 
volume for four representative values fh u = 10 -2 , 10 -3 , 10 -5 , and 10 -6 eV. In this analysis, 
we use the same values for the other parameters as in the analysis of Figure 2. For larger 
values of light neutrino mass {rh v > 1CT 3 eV), i.e. larger Yukawa coupling y u , the N decay 
term of Eq. (3.2) is strong enough to dominate other interaction terms of N before the 
annihilation effect of N becomes weaker than the dilution effect due to the expansion of 
the Universe. Therefore, the decay effect keeps N in thermal equilibrium for a longer 
time compared with the case that N is stable, and N can continuously remain in thermal 
bath before Ni is decoupled from thermal bath as can be seen from the top two panels 
of Figure 4. On the other hand, for smaller light neutrino masses (m u < 10~ 3 eV), the 
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Figure 3: Contour plots for the relic abundance of the right-handed sneutrino dark matter N% in 
the ttin ~ ra^, plane when the decay effect of iV is excluded. Each panel shows the cases = 
300 and 450 GeV, respectively. We fix the Z' gauge boson mass as Mz> = 1200 GeV. 

annihilation effect of N becomes weaker than the dilution effect before the ./V decay effect 
dominates other reactions. As a result, a retarded behavior appears in the number density 
evolution of N\ and N as can be seen from the bottom two panels of Figure 4. In addition, 
since the decay modes of N are governed by the Dirac neutrino Yukawa coupling y v , i.e. 
the effective neutrino mass fh u , the result in Figure 4 converges to the one in Figure 2 for 
smaller m v . 

3.2 Dependence on 

In this and the following subsections, we will obtain the thermal relic abundance of right- 
handed sneutrino O^/i 2 given mass parameters including m^,mN, fh u , m^,, Mz> and also 
the U(l)' gauge coupling g' by solving the coupled Boltzmann equations (3.1) and (3.2) 
including the decay term. In order to see the dependence of the relic abundance on each 
parameter, we search for the allowed regions, which satisfy the observed recent DM relic 
density limit [17], in the mj^-m^, plane fixing the other parameters. 

In Figure 5, the rriN-m^, parameter space is explored for the six cases of the right- 
handed sneutrino dark matter mass: = 150, 300, 450, 600, 750 and 900 GeV. Our 
numerical analysis is performed in the case of Mz> = 1200 GeV and rh v = 10~ 3 eV. The 
region between two thick dashed lines represents points where Q^h 2 is consistent with 
the recent WMAP result on the DM relic density [17]. As one can see in Figure 5, the 
right-handed sneutrino mass less than 900 GeV, < 900 GeV, is allowed by the recent 
WMAP result. 
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Figure 4: The actual number of N\ and N per comoving volume. The panels correspond to 
rh v = 10~ 2 , 10 -3 , 10~ 5 and 10~ 6 eV respectively from left to right and top to bottom. Blue solid 
and red dashed lines show =n^/s and Yjv = n^/s. The other parameters are the same as in 
the analysis of Figure 2. 

3.3 Dependence on m v 

The right-handed sneutrino DM particle N\ can remain in thermal equilibrium through the 
interactions of the right-handed neutrino N for which the inclusion of the decay effect is 
crucial as shown already. In this subsection, we study this effect in more detail by changing 
the neutrino Yukawa coupling y u , that is, the effective neutrino mass rh u . In Figure 6, the 
thermal relic density Q^h 2 is shown for the two representative cases m u = 10 _1 and 10~ 5 
eV in the ruN-m^, plane to be compared with Figure 5. We use the numerical values 
mft = 300 GeV and M z > = 1200 GeV for this analysis. 

For smaller effective neutrino mass, the decay rate of the right-handed neutrino is 
weaker, and consequently the right-handed sneutrino DM particle is decoupled earlier from 
thermal equilibrium. As a result, the relic abundance O^/i 2 increases as m u decreases. 
This implies that smaller Z' mass is required for smaller rh v . One can see the tendency 
from the top-middle panel of Figure 5 and Figure 6. 

3.4 Dependence on Mz> 

The right-handed sneutrino DM N\ is kept in thermal equilibrium through the interactions 
of the right-handed neutrino N which also depends on the Z' mass as the right-handed 
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Figure 5: Contour plots for the relic abundance of the right-handed sneutrino dark matter Ni 
in the mjv — rng, plane. The panels show the cases = 150, 300, 450, 600, 750 and 900 

GeV respectively from left to right and top to bottom. We fix the other parameters as follows: 
Mz' = 1200 GeV and fh v — 10~ 3 eV. In each panel, the region between two thick dashed lines is 
preferred by the recent result on the DM relic density. 



neutrino annihilates to the SM fermions trough the s-channel Z' mediation. In Figure 7, 
we therefore present the thermal relic abundance Q^/i 2 in the m^-m^, parameter space 
to check the dependence of the relic abundance Q^h 2 on the Z' mass Mz>- In the figure, 
the left-panel corresponds to Mz> = 2000 GeV and the right-panel to Mz> = 4000 GeV. 
In this analysis, we use the numerical values fn^ = 300 GeV and fh v = 10 -3 eV. The 
parameter space between two thick dashed lines is also the region which is allowed by the 
recent WMAP observational result on the DM relic density [17]. 

The decoupling of the right-handed sneutrino DM particle Ni from thermal equilibrium 
is determined by the N\N\ annihilation to right-handed neutrinos through the t-channel 
Z' exchange if the right-handed neutrino N remains in thermal equilibrium. Therefore, the 
relic density of the right-handed sneutrino DM barely depends on the U(l)' gauge boson 
mass Mz 1 if the right-handed neutrino remains in thermal equilibrium by some other 
interactions when the DM particle N\ is decoupled. This is the case with large values of 
the effective neutrino mass, fh u > 10 -3 eV, for which the decays and inverse decays of the 
right-handed neutrino become dominant interactions before the annihilation of N freezes 
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Figure 6: Contour plots for the relic abundance of right-handed sneutrino dark matter Aq in the 
to jv — 777^, plane. Left and right panels are respectively corresponding to 777^ = 10 _1 and 10 -5 eV. 
The parameters and mz< are fixed as = 300 GeV and Mz> = 1200 GeV. In each plane, 
the region between two thick dashed lines is allowed by the recent DM relic density result. 

out. Thus, the right-handed neutrino remains continuously in thermal equilibrium during 
the freeze-out of the DM Aq as can be seen clearly from the top-left panel of Figure 4 with 
fh u = 10~ 2 eV. Consequently, one can expect that the relic abundance of Aq has almost 
no dependence on the Z' gauge boson mass. However, for smaller values of the effective 
neutrino mass, m u < 10~ 3 eV, the annihilation of N may freeze out before the decay 
modes of N become important, and hence the decoupling effect of the AW annihilation 
is non- negligible. As Mz> increases, the AW annihilation cross section becomes smaller, 
and thus the AW annihilation decouples earlier. Because of this, the relic density of the 
right-handed sneutrino DM becomes larger. Comparing the top-middle panel of Figure 5 
and the two panels of Figure 7 with Mz 1 = 1.2,2 and 4 TeV, respectively, one can find 
the fact that the thermal relic abundance increases as the Z' gauge boson mass 

increases, and thus the right dark matter density is obtained for smaller Z' mass. 

3.5 Dependence on g' 

In all analysis of this paper, we use g' = a/5 / 3 #2 tan 0w ~ 0.46 as the reference value for 
U (1)' gauge coupling [7]. However, in this section, we vary the U(l)' gauge coupling to show 
the dependence of the relic abundance of Aq on this gauge coupling. In the left two panels 
of Figure 8, we present the thermal relic density of Aq for two cases g' = 0.3 and g' = 0.25 
with mfr = 300 GeV; in the right panel, for the case g' = 0.2 with = 150 GeV. In all 
the cases, we use the same numerical values for the other parameters: Mz> = 1200 GeV 
and m v = 10~ 3 eV. As in the previous cases, the regions between two thick dashed lines 
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Figure 7: Contour plots for the relic abundance of the right-handed sneutrino dark matter Ni 
in the mjv — rug, plane for the cases Mz> = 2000 and 4000 GeV. The parameters and m y 
are fixed as = 300 GeV and fh v = 10~ 3 eV. The regions between two thick dashed lines are 
preferred by the recent result on the DM relic density. 
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Figure 8: Contour plots for the relic abundance of the right-handed sneutrino dark matter N\ in 
the rriN — m%, plane. The left two panels respectively show g' = 0.3 and <?' = 0.25 with — 300 
GeV, and the right panel shows g' — 0.2 with = 150 GeV. We fix the other parameters as 
follows: M z < = 1200 GeV and m v = 10^ 3 eV. The regions between two thick dashed lines are 
allowed by the recent DM relic density observation. 

are consistent with the recent DM relic density observation [17]. 

Obviously, for smaller g' , the interactions of N\ and N become weaker, and thus the 
right-handed sneutrino DM is decoupled earlier from thermal bath. Consequently, the 
thermal relic density of DM increases, which can be seen from the top-middle (top-left) 
panel of Figure 5 and the left two panels (right panel) of Figure 8. As can be seen from the 
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left two panels of Figure 8, the right-handed sneutrino is overproduced for g' < 0.25 if N\ 
is heavier than 300 GeV. When N\ is light enough, fn^ ~ 150 GeV, the allowed parameter 
space marginally exists even for g' = 0.2. Such a low limit is almost independent of the 
effective neutrino mass if m u > 10 -3 eV. 

As a result, an interesting lower bound on the U(l)' gauge coupling g 1 > 0.2 can be put 
by the current dark matter relic density limit [17] in the U (l) x model with the right-handed 
sneutrino dark matter. 

4. LHC signatures 

For the smoking gun signal of this model, we look for the productions and decays of the 
new particles in the model. These include the extra gauge boson Z' , its superpartner Z' ', 
the right-handed neutrino N and the LSP dark matter N± being a superpartner of N. 
The standard search for an extra gauge boson Z' is through the observation of high mass 
dilepton resonances from Z' — > l + l~ [9]. When the extra U(l)' is associated with the 
seesaw mechanism as in this model, it is also important to look for the mode Z' — > NN 
[10, 11, 12] to test the Majorana nature of neutrinos and the seesaw mechanism. In this 
section, we will discuss the associated phenomenology of Z' ', which is required to be lighter 
than Z' as was discussed in the previous section, in parallel with the Z' phenomenology. 

4.1 Production and decay of Z' 
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Figure 9: The Z' production cross-section times the branching ratios for the SM channels including 
right-handed neutrinos. Two panels are for V§ = 7 and 14 TeV, respectively, with tan/3 = 10. 



Single Z' produced at the LHC like the Standard Model Z can be observed through its 
decay modes. All the possible decay channels are shown in Table 1 for the supersymmetric 
U(l) x model. Our main focus is the production of the right-handed neutrinos through 
the process pp — > Z' — > NN which can be compared with the standard Z' discovery 
channel pp — >■ Z' — )■ I P. In this paper, we assume that the Z 1 decays to sfermions are not 
allowed kinematically. In Figure 9, we show the production cross-section for each SM decay 
channel of Z' including right-handed neutrinos in the decoupling limit of non-SM Higgs 
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bosons. The mass of the right-handed neutrino was taken to be 300 GeV for the set of 
plots. The renormalization/factorization scale was chosen to be V~§ and CTEQ6L [18] was 
taken as PDF for the cross-section calculation. Apart from the dilepton and NN channels, 
the Higgs modes could also be interesting. For the quark, lepton and right-handed neutrino 
modes in Figure 9, we have shown only one flavor contributions. One can see that the down 
type quark mode has larger cross-section than the up type quark as well as the leptonic 
modes, which is consistent with the charges given in Eq. (2.1). By the same reason, the NN 
production cross-section is the largest, which enhances the Majorana neutrino signatures 
(see Subsection 4.3) compared to other U(l)' models. Concerning the Higgs contribution, 
we have considered only the SM Higgs which leads to the Zh mode. This has a very low 
cross-section for large tan /3 as the coupling is proportional to cos(q + /3) which is sin(2/3) 
in the decoupling limit. Thus, for low tan/3, the Higgs cross-section gets enhanced and 
comparable to the fermionic modes. Let us also remark that the production cross-section 
of the standard Z' discovery channel, Z 1 — > l^P, depends on tan /3 and on the number of 
new channels, e.g., with Higgs bosons or sfermions that can be kinematically allowed in 
our model. In the case of a U(l)' extension of the non-supersymmetric Standard Model, 
the model dependence comes only from the number of right-handed neutrino channels. 
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Figure 10: Same as in Figure 9 but with non-decoupling heavy Higgs bosons. 



Figure 10 describes the case with all the possible Higgs decay modes, which includes 
Z' decays to HZ, hA, HA and H^H T . For this study, we have taken the CP-odd Higgs 
to be non-decoupled , i.e. tua = 162.2 GeV, which makes the charged Higgs mass 180 
GeV and the heavy neutral Higgs 164 GeV. As can be seen in Figure 10, the production 
cross-sections are not promising at least for the case of center of mass energy of 7 TeV. In 
the 14 TeV case, H±E* and hA can approach the production cross-section ~ 100 fb at 
tan/3 = 10. The other Higgs channels, HA, hZ and HZ, do not have enough effective-cross 
section at 14 TeV. At tan/3 = 10, sin(a + /3) ~ 1, which makes the h-A-Z' and Z'-^-H^ 
couplings almost same (see Table 1) which is also reflected in Figure 10. The ratio of the 

production cross-section between hA and HA is proportional to ( ~(~q^ J and the same 
is the case for ZH and Zh, as is expected from the nature of their couplings given in Table 1. 
It is also to be noticed that Higgs pair production channels have more production rates 
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compared to associated gauge boson channels due to the momentum dependent couplings 
of the formers. Given the suppressed Higgs channels, the production cross-sections for the 
NN and 11 modes remain almost same as in the decoupling limit. 

4.2 Production and decay of Z' 

Unlike Z' whose mass is strongly constrained by the precision electroweak data, its super- 
partner Z' can be light, as required in the model of the right-handed sneutrino dark matter. 
Thus, it can lead to interesting phenomenology related to the seesaw mechanism even 
though Z' turns out to be too heavy to be produced at the LHC.Let us first try to look for 
the direct pair production of Z' at the LHC. In Figure 11, we estimate the pair production 
rate, i.e. pp —> Z'Z' , with the variation of m^, at the LHC for the center of mass energy 
of 7 TeV and 14 TeV. The renormalization/factorization scale was chosen to be y/§ and 
CTEQ6L [18] was taken as PDF for the cross-section calculation as was in the previous 
case. The rates appear to be very low for both 7 and 14 TeV cases, due to the fact that only 
t-channel electroweak diagrams contribute to the process. This results in the production 
cross-section ~ fb even for the light Z' . 

10 
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Figure 11: Variation of Z' Z' production cross-section. 

However, let us note that the LHC being a machine with huge gluon flux, the strongly 
interacting supersymmetric particles namely squarks and gluino can be copiously produced. 
Cascade decays of such supersymmetric colored particles could then be a good source of 
superpartners of the electroweak particles. In particular, squarks decay through the direct 
electroweak couplings to quarks and neutralinos or charginos, and gluinos decay through 
the two-body decay to squarks and quarks or through three-body decay to quark pair and 
charginos or neutralinos. An interesting application of such supersymmetric cascade decays 
to the Higgs production has been studied [19] in the context of MSSM. 
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Figure 12: Variation of strong production cross-section. 



Figure 12 describes the variation of strong production rate with the common mass 
M str0 ng of squarks and gluinos for the center of mass energy of 7 TeV and 14 TeV including 
the contributions from only first two generations of squarks. For M str0 ng = 1 TeV, the 
cross-sections are 43 fb and 1.3 pb for the 7 and 14 TeV center of mass energy, respectively. 
The inclusion of third generation will of course enhance the production rate. This can be 
compared to the Z' — >■ NN production cross-sections for Mz> = 1 TeV: 70 fb and 350 fb 
for the 7 and 14 TeV for one generation of right handed neutrino, respectively. Inclusion 
of three generation will increase the production cross-section by factor 3. The numbers of 
Z' and Z' events will sensitively depend on the masses of squarks/gluinos and Z' . The 
squarks and gluinos can now decay to Z'. The branching fraction to Z' depends on the mass 
parameters and on the coupling with left and right-handed up and down type quarks as 
shown in Eq. (2.1). Figure 13 describes the q — > qZ' branching fraction with the variation 
of the squark mass for rug, = 300 GeV, 500 GeV and 700 GeV. From the figure, we can 
see that the du — >■ dZ' has larger branching fraction compared to other quark modes as is 
evident from Eq. (2.1). The charge corresponding to dji type quark is 3/2\/l0 which is the 
largest among quarks. Note that the branching ratio of q — > qZ' reaches about 70% unless 
the kinematic suppression is applied. This will lead to a sizable number of Z' produced 
from the cascade decays of squarks and gluinos. In the case of Z' being the next lightest 
supersymmetric particle (NLSP), all of the pair produced strong particles will end up with 
a pair of Z', i.e. the effective branching fraction is 100 %. 

The decay of Z' can give rise to various final states depending on the possible decay 
modes that are open: that is, Z 1 — > NN ,11, HH and S12S12. Here H and H denote any 
type of Higgs bosons and Higgsinos shown in Table 1. The U(l)' Higgs bosons, Sip, or 
Higgsinos, Sip, are expected to be as heavy as Z' and thus heavier than Z'. Assuming 
the Higgsinos heavier than Z' , we will concentrate on the first two decay channels in this 
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Figure 13: Variation of Br(g — > qZ') with M strong for m*, = 300, 500, and 700 GeV respectively. 



paper. Note again that the Z' — > NN mode is the dominant mode given the U(l) x charge 
assignment (2.1) and becomes the unique one if Z' is the NLSP and N is the LSP, N\. In 
the next subsection, we will focus on the mode Z' — > NN\. 



4.3 Signatures of the seesaw and displaced Higgses 

As discussed in the previous subsections, the extra U(l) x gauge boson Z' and its super- 
partner Z' can be copiously produced and decay to right-handed neutrinos at the LHC. 
Thus, we can search for the signatures of the pair produced heavy Majorana neutrinos 
through two channels: 

pp^Z'^NN, (4.1) 
pp -> Z'Z' -> iViViViA>i . 

Now, the right-handed neutrino, N, can decay to the following final states: 

N -»• IW, vZ, vh, vH, uA, IH + (4.2) 

as shown in Subsection 2.2. 

Figure 14 describes the mass variation of the decay branching fraction of the right- 
handed neutrino and the partial decay widths for various decay modes for the effective 
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Figure 14: Branching fractions and partial decay widths of the right-handed neutrino, N, in the 
non-decoupling limit of heavy Higgs bosons with = 162.2 GeV as a function of for the 
effective neutrino mass, m v = 1 meV and tan j3 = 10. 



neutrino mass rh v = 1 meV. From Figure 14, it is clear that the right-handed neutrino 
decays to gauge bosons and lepton modes have more decay branching fraction ~ 35 — 45% 
for the whole region of the right-handed neutrino mass, uin- Apart from the production 
of the gauge bosons from the decay of the right-handed neutrino, Higgs modes are also 
possible due to the direct coupling in the superpotential Eq. (2.2) which is proportional to 
the small Yukawa coupling, y v . The mode N — >■ hv has larger branching fraction compared 
to the ./V — > Hv mode, which is also expected from Eq. (2.18). The right-handed neutrino 
decays to the pseudo-scalar Higgs mode, i.e. iV — > Au, and also the charged Higgs modes, 
i.e iV — > H^P, have very low branching fraction because of the fact that the decay widths 
are proportional to cos 2 (3 ~ 10 -2 . For the choice of smaller tan/3, these modes also could 
be interesting. 

In the model under discussion, the Majorana nature of right-handed neutrino, N, can 
be probed either of the two channels (4.1) leading to the same-sign dilepton (SSD) final 
states: 

PP -> Z' (Z'Z') -> iH+W^W* (+ jfr) . (4.3) 

Depending on the decay of W^, the final state can have 3/ or 4/, or SSD + 4j (+ r/ T ). The 
missing energy contribution comes from Z' decaying to the LSP, N\. 

It is also interesting to look for the Higgs signal from the channel N — >• hv. When 
the Yukawa coupling, y u , or the effective neutrino mass, fh u , is small enough, the Higgs 
thus produced will be displaced and its main decay to bb can be observed [13]. If the other 
right-handed neutrino, N, decays to V^W^ or vZ, then the final states can have one or 
two charged leptons to tag along with the bb: 

pp^ Z' , Z'Z' -> h l ± W T /Z+ p/ T . (4.4) 

To get some reference values for the production cross-section of this signal, let us take 
Br(N l^W^) x Br(N — > hv) ~ 5% from Figure 14 and there is a combinatorial 
factor 2 as one of the N has to decay to Higgs which makes Br(NN — > hvl^W^) ~ 
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10%. At the 7 TeV LHC, we have a(pp -> Z' ->• iVJV) ~ 0.07 pb for M z / = 1 TeV, and 
a(pp -> Z'Z' -> iVJV) = 43 fb for M stro „ 9 = 1 TeV assuming the Z' NLSP as was discussed 
in the previous subsections. This leads to the production cross-section of the process (4.4): 
a(hl^W Zfl ) = 21 fb and 4.3 fb from the Z' and Z' channel, respectively. Thus, there is a 
chance to find the Higgs signal at the 7 TeV LHC if, in particular, the associated displaced 
vertex is large enough to kill the backgrounds. The corresponding figures at the 14 TeV 
LHC are aihl^W^^) = 105 fb and 130 fb from the Z 1 and Z 1 channels, respectively. 

Figure 15 describes the variation of the decay length of right-handed neutrino, N, with 
the right-handed neutrino mass, rn/v, for two effective masses rh v = 1 and 0.01 meV. One 
can see that the decay length becomes favorably large for rh u < 1 meV making clean the iV 
decay signals. There is a large parameter space for such a effective neutrino mass allowing 
the right-handed sneutrino dark matter as shown in Figure 5. 
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Figure 15: Variation of the decay length of the right handed neutrino for m v = 1 and 0.01 meV, 
respectively. 



For small tan/3, the right-handed neutrino, N, can have a large branching fraction to 
the charged Higgs and charged lepton, if kinematically allowed. Thus, we can have SSD 
plus charged Higgs in the final state: 

pp -> Z' (Z'Z 1 ) -> H^W^FF (+ pr) (4.5) 

in which the charged Higgs decays to tv t or tb. These lead to the displaced multi-jet (r-jet 
or 6-jet) and multi-lepton final states. 

When Z' is not the NLSP, Z' can also contribute to 11* if kinematically allowed. In 
this situation, we can have the following final states: 

( W ± Pl ± P+p / T 
Z°l ± P+ fa 



pp 



Z'Z' 



(4.6) 
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which may involve one displaced and one prompt vertex. 

The final states discussed above can be studied as smoking gun signals for this model 
at the LHC. Among those channels, the channels associated with the multi-lepton final 
states suppress the SM backgrounds effectively [20]. Again the displaced decay of the 
right-handed neutrino actually can remove the SM background completely. Nevertheless, 
still for the estimation of signal significance, one needs to consider the following back- 
grounds: ti, tiZ, tih, W + n— jets, and Z + n— jets as well as other supersymmetric final 
states. The detailed simulation for the above mentioned final states and the corresponding 
significance calculation deserve further investigation [21]. There is another aspect of this 
model that could be interesting through the mixing in the neutrino (2.15) and sneutrino 
(2.17) sector. Specially, in the context of the NLSP decaying to the LSP (Aq), these can 
lead to remarkable features in the final state depending on the nature of the NLSP [21]. 

5. Conclusion 

We considered the possibility of a right-handed sneutrino as the LSP dark matter in the 
supersymmetric Standard Model extended to include an extra U(l)' gauge symmetry real- 
izing the seesaw mechanism with three right-handed neutrinos. In a supersymmetric seesaw 
model, a complex right-handed sneutrino gets split into two real mass eigenstates due to the 
soft supersymmetry breaking Majorana mass term. While the lightest real right-handed 
sneutrino, Aq as the LSP, cannot annihilate through the U(l)' gauge boson Z' exchange at 
the s-channel, its annihilation to lighter right-handed neutrinos, N\N\ — > NN, mediated 
by the U(l)' gaugino Z' at the t-channel is shown to be effective in generating the right dark 
matter relic density. In this process, the decay and inverse decay of N play an important 
role in maintaining N and Aq longer in thermal equilibrium and thus reducing the dark 
matter density. This behavior was shown in Figure 4 by solving the Boltzmann equations 
with varying the decay rate quantified by the effective neutrino mass, rh u . The resulting 
dark matter density is computed in Figures 5, 6, 7, 8 showing the favorable parameter 
space of the model. 

The signatures of the model can be probed through the production of Z' as well as Z' 
and the corresponding decays. In particular, the production of the right-handed neutrino 
is of great interest because of its Majorana nature and possible displaced vertices. When 
Z' is light enough as is needed for the right-handed sneutrino dark matter, it can be pair- 
produced copiously through the cascade decays of squarks/gluinos as estimated in Figures 
12,13 leading to a large number of events for the process of pp — > Z'Z' — > NN+ pr- This 
can be compared with the production cross-section of the usual process pp — )■ Z' — > NN as 
shown in Figures 9, 10 including all the other final states except the sfermions. The above 
channels provide the golden search for the seesaw mechanism and the Majorana nature 
of neutrinos through the same-sign dilepton final states. We also point out a remarkable 
feature of the Higgs production from the right-handed neutrino decay. The displaced bb 
along with displaced tagged leptons will be a clean signature in probing the light neutral 
Higgs. The non-decoupled heavy Higgs bosons can also be probed in a similar way. In 
this article, we have given the number of the production rates and the effective branching 
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fractions which gives a hint that these signal topologies can be probed with early data of 
the LHC for the 7 TeV case. The detailed collider simulation to calculate the acceptance 
under the basics and the hard cuts will be reported in a separate work [21]. 
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